Rectification of self-propelled particles by symmetric barriers 
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The motion of self-propelled particles can be rectified by asymmetric or ratchet-like periodic pat- 
terns in space. Here we show that a non-zero average drift can already be induced in a periodic 
potential with symmetric barriers when the self-propulsion velocity is also symmetric and period- 
ically modulated but phase-shifted against the potential. In the adiabatic limit of slow rotational 
diffusion we determine the mean drift analytically and discuss the influence of temperature. In the 
presence of asymmetric barriers modulating the self-propulsion can largely enhance the mean drift 
or even reverse it. 
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Ratchets and rectification of passive Brownian par- 
ticles or the idea of generating "noisy transport far 
from equilibrium" have been a central topic of statistical 
physics for the last 20 years with far reaching impact on 
biology and in applications [H-Q. In particular, a detailed 
analysis under what conditions and symmetries a ratchet 
allows directed transport revealed its basic concepts [l| . 

Recent years have seen a tremendously growing in- 
terest in understanding self-propulsion in nature and of 
artificial microswimmers 0, Exploring the conse- 
quences of self-propulsion currently develops into another 
paradigm of noncquilibrium statistical physics (see, for 
example, 0-0] )• Also here, rectification of self-propelled 
particles by static asymmetric barriers Tol- 13j has re- 
cently attracted much attention due to prospective ap- 
plications of this easily realizable phenomenon in nano- 
and micro-scale devices. Experimental studies with E. 
coli supported by theoretical work 0, Hi -17 1 show the 
pivotal role of self-propulsion for rectifying cell motion 
in an array of asymmetric funnels IIOL 111! or for driving 
a nano-sized ratchet-shaped wheel [12i.ll3l|. Our goal here 
is to put forward a novel concept of the rectification of 
active particles, when one is able to artificially modulate 
their self-propulsion velocity. 

Janus particles and their anisotropic surface proper- 
ties are ideal for generating active motion, for example, 
through self diffusiophoresis in catalytic reaction prod- 
ucts II 811 or self thermophoresis under illumination with 
light HH. InRefs. [H[21 



lated v(z) on a driven pulsating ratchet of a passive parti- 
cle. We then demonstrate analytically that in a periodic 
potential with symmetric barriers active particle motion 
can be rectified with the help of a modulated propul- 
sion velocity. Furthermore, in potentials with asymmet- 
ric barriers, a modulated v(z) can largely enhance the 
mean particle drift and even reverse it. 

We consider an active Brownian particle that moves 
in a spatially patterned environment with a one- 
dimensional periodic structure characterized by a po- 
tential U(z) with periodicity L. In addition, we as- 
sume the propulsion velocity of the particle v(z) to be 
also periodically modulated in space along the same di- 
rection as the periodic pattern. More specifically, we 
set v(z) = vq[1 + w(z)], where vq is the mean self- 
propulsion velocity. The modulation function w(z) has 
the same periodicity L and satisfies J" + w(z) dz = 
and min[l + w(z)} > 0. 

In the overdamped limit, the position r = (x,y,z) 
of the active particle and its intrinsic direction of self- 
propulsion given by the unit vector p obey the Langevin 
equations 
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the Janus particles are placed 
into a critical fluid mixture which allows to fine-tune their 
self-propulsion velocity v by the intensity of the applied 
light field. In particular, with scanning optical tweezers 
one can then create a periodic modulation of v in space 
in addition to the presence of a periodically corrugated 
channel or a set of artificially created barriers [i(J ■ 

Motivated by such an experimental system we show 
in this letter that an arbitray modulation of the self- 
propulsion velocity v(z) alone does not generate a mean 
particle drift. We map the overdamped active Brownian 
particle moving in a static potential and with a modu- where we introduced the L-periodic and time-dependent 



r = -fj,d z U(z)e z + v(z)p + £(t), p = rj(t) x p. (1) 

Here, \x is the translational mobility and the unit vec- 
tor e z points along the z-axis. and r/(i) rep- 
resent translational and rotational noise, respectively, 
with zero mean. Their variances are linked to the 
translational (fiksT) and the rotational (D r ) diffusivi- 
ties via the fluctuation-dissipation relations (£(i)£(f')) = 
2fj,k B TlS(t - t') and (r](t)r)(t')) = 2D r l6(t - t'). 

According to Eqs. ([1}, the motion in the x,y plane is 
decoupled from the motion along the z-axis. Introducing 
q = cos 8 = p ■ e z , we can rewrite the equation for the 
z-coordinate of the particle in the form 



dU cS {z,q{t)) 

v a q{t) h£(i), 

oz 



(2) 
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effective potential U e s(z,q(t)) 

U eS (z,q(t)) = fiU(z) + q(t)W(z), 



(3) 



with W{z) = —vq J w(z) dz. One can show that the 
stochastic process q(t) satisfies the Stratonovich equation 



q = -D r q + ^l-q 2 ri{t), 



(4) 



where rj(t) is a Gaussian white noise with (r)(t)r](t')) = 
2D r 5(t — t'). The stationary distribution R(q), and the 
time correlation function (q(t)q(t')) are given by 



R(q) = 1/2, (q(t)q(t')) = (l/3)e' 



-2D r \t-t'\ 



(5) 



respectively [231 ] . 

Equation ([2]) shows that an active Brownian particle 
moving in a patterned environment with a spatially mod- 
ulated propulsion velocity is equivalent to a passive Brow- 
nian particle moving in a randomly pulsating potential 
U e ff(z, q(t)) and driven by a random force Voq(t). The lat- 
ter is widely known 'driven pulsating ratchet" 0,[H. 
Here we have a special type of stochastic pulsation and 
a stochastic driving force oc Voq(t). 

Now we investigate under which conditions the active 
particle system generates a non-zero mean drift veloc- 
ity (i). First, using similar arguments as in Ref. pij . we 
show that in the absence of the external potential, U(z) = 
0, the average drift of the particle vanishes exactly for 
an arbitrary spatially modulated propulsion velocity v(z) 
meaning (i) = 0. We introduce an auxiliary stochas- 
tic process Z{t) = z{—t) with Z = —dz(—t)/d(—t) and, 
consequently, obtain (Z) = —{z). The time evolution 
equation for Z(t) follows from Eq. ([2]), 



Z = - Vo q(-t)-w(Z)q(-t)-Z(-t). 



(0) 



We recall that — £(— t) is again a Gaussian white noise. 
Also the statistical properties of —q(—t) = —e z -p(—t) are 
identical to those of q(t) (parity-time invariance since 
p(t) simply performs a random walk on the unit sphere. 
Consequently, we can substitute in Eq. (J6j) —£(—t) by 
and —q(—t) by q{t) and obtain (Z) = (i) = — (i) = 0. 

Second, in a symmetric potential U(z) the active- 
particle current is zero for constant self-propulsion ve- 
locity. However, we now show that already a symmetric 
velocity modulation w(z) can be suffcient to induce a 
non-zero mean current. Again, we introduce an auxiliary 
variable Y{t) = —z(t) + zo, with some arbitrary displace- 
ment zq and obtain from Eq. ([2|), 



Y = -v Q q(t)-w(-Y(t)+z Q )q(t) 
+ fjU'(-Y{t) + ZQ)-^t). 



(7) 



We replace the stochastic processes — q(t) and — in 
Eq. (0 by q(t) and using their time- reversal symme- 
try. We then ensemble average Y and conclude that the 
mean current is always zero, (Y) = (i) = — (z) = 0, when 




UJ-q,z) 



FIG. 1: (Color online) Schematic representation of the rec- 
tification. The active particle experiences different effective 
potentials Eq. ([3]) (shown by shaded area), when moving in 
opposite directions. 



a unique displacement zq exists, for which both potential 
U(z) and velocity modulation w(z) become symmetric 
about z = 0. This means U(—z + zq) = U(z + zq) and 
w(— z + zq) = w(z + Zq). So, only when the left-right 
symmetry along the z direction is broken, for example, 
by a non-zero phase shift between U{z) and w(z), can 
one obtain a non-zero mean current. Another possibility 
would be to introduce a ratchet-like modulation of w(z). 
We concentrate here on the first case. 

Symmetric potential and velocity modulation: We now 
investigate in more detail the most striking manifesta- 
tions of the last result that one can rectify active swim- 
mers with a symmetrically modulated self-propulsion ve- 
locity in a static and symmetric potential. For the sake 
of argument, we set v(z) = vq[1 + a sin (kz + 5)] and 
U(z) = Uq sinfcz, where we introduced the wave number 
k = 2n/L. From the above symmetry arguments we ex- 
pect a non-zero current (z) ^ for 6 ^ 7m, n = 0, ±1, . . . 

In order to unveil the bare physical mechanism of the 
rectification we will carry out the following analysis for 
zero temperature, i.e. £(t) = in Eq. ([2]). In addition, 
we will assume that the stochastic process q(t) changes 
adiabatically slow in time. In case of active Brownian 
particles this corresponds to small rotational diffusivities. 
More precisely, the rotational decorrelation time oc D^ 1 
should be much larger than the time L/vq needed for the 
active particle to swim the distance L of one period of 
the potential. 

In Fig Q] we illustrate the basic mechanism for the rec- 
tification, where we plot the effective potential U e s for 
particles moving in positive (+q > 0) and negative (—q) 
direction. Clearly, the particles moving to the right ex- 
perience steeper potential rises or larger opposing forces, 
which can even trap the particle in a mimimum if the 
force is too large. Whereas in the opposite direction, 
the potential here is nearly constant. So, the mean drift 
velocity (z) will be negative. Only for phase differences 
5 = 0, 7T, . . ., the shapes of U c s{+q) and U c s(—q) are iden- 
tical but phase shifted against each other and (z) = 0. 

With the potential U(z) and the velocity modulation 
w(z) both sinusoidal, the noiseless Eq. can be conve- 
niently written as the well-known Adler equation pjj 

z = v q + A(q) sin [kz + (f){q)] , (8) 

with amplitude A{q) = \J (qavo cos S) 2 + (qavo sin 6 — fc/i[/o) 2 
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and phase shift tan</> = (qavosinS — kfillo)/qavocos8 
considered as a constant in the adiabatic limit. 

The long-time solution of Eq. (JSJ is either a pair of two 
fixed points in z £ [0,L] for \q\vQ < A(q) meaning the 
active particle with z = is trappend in a minimum of 
Uctf — voqz. Or the time-dependent solution for \q\vo > 
A{q) increases (q > 0) or decreases (q < 0) in time. In 
fact, the solution for |g|-uo > A(q) can be decomposed into 
a periodic component and a dc component representing 
the average slope (z)(q) (time averaged z at fixed q), 
which is given by 

(z) (q) = qv y/1 - (A(q)/qv ) 2 . (9) 

The average zero temperature (zero-T) current (i) is 
found by averaging the slope (z) (q) from Eq. (|9]) over 
all possible q with the stationary distribution function 
R(q) = 1/2 from Eq. ©. i.e. as (i) = J^(z)(q) R(q) dq. 

For arbitrary phase shift 6 and amplitude a of the ve- 
locity modulation, (z) can be computed analytically but 
is given by a rather lengthy algebraic expression. For the 
sake of clarity, we only mention the zero-T current in two 
special cases: (i) when a = 1 and 5 = tt/2 (meaning the 
self-propulsion is zero at specific points in space), and 
(ii) when vq — > oo. The corresponding zero-T currents 
(i)W'( M ) are given by 

w = VMo {2vq _ W 3/ 2 {i>W = fcM/osin£ (I0) 
D^o VI — a 2 

From Eqs. (JTUJ) we conclude that the current remains fi- 
nite in the limit vq — > oo for | a |< 1, and that for 
a = 1, the current increases unlimited with vq accord- 
ing to the asymptotic law (i) ~ y/kfiUovo- Furthermore, 
by inspecting Eqs. ((HJ and ([9]), we find that the optimal 
rectification regime is reached for S = tt/2 and a = 1. 

The system studied so far has striking similarities with 
the classical Seebeck ratchet, where a passive particle 
moves in a static periodic potential under a periodic tem- 
perature profile [l[. Typically, both fields are also sym- 
metric and have the same periodicity. Only when the spa- 
tial symmetry of the system is broken, i.e., when there is 
a non-zero phase shift between potential and temperature 
profile, does one observe a mean particle velocity exactly 
as in the case discussed here. The profound difference is 
the term with the self-propulsion velocity vo in Eq. 
Even at zero temperature and in the limit D r — > 0, where 
q only changes adiabatically, can one observe a non-zero 
current of active particles. 

For any non-zero temperature, the rectification cur- 
rent (i) of active Brownian particles can be found in 
the adiabatic approximation by using standard results 
for passive Brownian particles 0]. Thus, for any fixed 
value of q = q(t), the current J(q) is obtained as the 
current of a passive particle that moves in the effective 
potential Eq. Q and is driven by the constant effective 
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FIG. 2: (Color online) Average current (J) in nondimensional 
units for U (z) = sin 2ttz and w(z) — a sin (2ttz + S) with D r = 
20. Solid red and dashed lines show the zero-T and finite-T 
adiabatic limits of (J), respectively. Symbols correspond to 
numerical solutions of the dimensionless Eqs. {lj. (a) (J) vs 
vq for different a and with 5 — n/2 and T = 1. Inset shows 
the 8 dependence at a = 1 and vq = 60. (b) (J) vs vq for 
different temperatures T at a = 0.5 and S = tt/2. Symbols 
belong to D r = 400 and T = 1. 

force v q. According to [l| we find the stationary condi- 
tional density p s {z \ q) and the current J{q) and give their 
lengthy formulas explicitely in the supplemental mate- 
rial. In what follows, we refer to J(q) as the finite-T adi- 
abatic current. The total rectification current (J) then 
becomes (J) = | J_ x J(q) dq. 

In order to illustrate our analytical results, we compare 
the finite-T and the zero-T average current (J) in the adi- 
abatic limit to numerical results by solving Eqs. ([TJ. We 
use reduced units, where we rescale position z by period 
L and time t by the characteristic time r = (/it/o/ T 2 ) -1 
which a passive particle needs to travel a distance L when 
it is driven by force Uq/L. We also refer all energies to Uq 
so that in particular the potential amplitude becomes 1. 
This introduces nondimensional parameters for which we 
use the same symbol; namely the reduced mean propul- 
sion velocity vqL/^lUq — > Vq, the reduced temperature 
ksT/Uo — > T, and the rotational diffusivity tD,, — > D r . 

In Fig. [2] we show several graphs of the average current 
(J). Solid red and dashed lines always correspond to 
the zero-T and finite-T currents in the adiabatic limit. 
Symbols are obtained from direct numerical solutions of 
the nondimensional Langcvin equations ([1]) with D r ~ 
20, implying that the adiabatic limit is expected for vo ^ 
20. 

In Fig.[2][a) both the finite-T and the zero-T currents 
increase monotonically with vo and saturate at the re- 
spective limiting value (z)^ given in Eq. (|T0| for any 
| a | < 1. The overall effect of temperature is to reduce 
the current since it disturbes the directed motion of the 
active particles and thereby decreases the difference be- 
tween particle currents in opposite directions. As pre- 
dicted by (i)W in Eq. (fTUl) . the zero-T current diverges 



4 



at a = 1 as vo — > oo. This implies that the current can be 
made as large as one desires by increasing vq. However, 
the current grows rather weakly since (J) = (i) ~ ^/vq 

The finite-T current increases monotonically with a. 
Moreover, (J) is a 27r-periodic function in 6 and symmet- 
ric about 5 = 7r/2 as the inset of Fig.^a) demonstrates. 
So the maximum current occurs at S = tt/2. 

Figure^b) clearly shows that the zero-T current is 
zero until vq reaches a value of ca. 2tt, the maximum 
force from the static potential U(z). The active particle 
needs a non-zero vq to move against the potential force. 
Temperature fluctuations help the particle to cross po- 
tential barriers and thereby activate a non-zero current. 
As already mentioned, for vq > 2tt temperature fluctua- 
tions smear out the difference between the active motion 
in opposite directions, and, as a consequence, the zero-T 
current is, generally, larger than the finite-T current. 

Interestingly, large temperatures such as T = 10 in 
Fig.[2£b) no longer initiate a non-zero current at small 
vq since they dominate the active-particle motion com- 
pletely. Furthermore, the largest rectification current is 
generally achieved in the adiabatic limit since rotational 
diffusion merely diminishes the current as demonstrated 
in Fig(5Jb) with numerical results at D r = 400. 

Strong current enhancement and reversal: Now, 
we demonstrate that a weak modulation of the self- 
propulsion velocity can have a pronounced effect on the 
average active current, which is induced by a static asym- 
metric or ratchet-like potential. The modulation either 
very strongly enhances the current or even reverses it. To 
this end we set U(z) = sin27rz+(l/4) sin47rz (see inset of 
Fig. [3} and use the same sinusoidal velocity modulation 
function w(z) as before. 

In Fig. [3] we plot (J) versus the amplitude of the ve- 
locity modulation a, which we obtained in the finite-T 
adiabatic limit for Vq = 10 which is just above the maxi- 
mum slope of the static potential, max(dU (z) / dz) = 9.3. 
We scale the current in units of the reference current 
(J)(o = 0) when the self-propulsion velocity is not mod- 
ulated. Different values of the phase shift parameter 5 of 
the velocity modulation function w(z) are used. 

The strongest effect of the velocity modulation w(z) 
on the average current is found for S = tt/2 and when 
a = — 1 meaning that the self-propulsion vanishes at 
specific points in the potential and doubles in between. 
In this case, the current is enhanced by more than an 
order of magnitude. Much more sensitive with respect 
to the modulation amplitude a is the sign of the average 
current. We find current reversal at any value of the 
shift parameter 6. Remarkably, for 5 — tt/2, the current 
is reversed by varying the propulsion velocity by as much 
as 5%, i.e. a = 0.05, as shown by the red line in Fig. [3] 

Controlling the motion of active or self-propelled Brow- 
nian particles is of immense interest both for technologial 
applications but also for a fundamental understanding. 
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FIG. 3: (Color online) Finite-T average current (J) in units of 
{J)o = {J)( a = 0) as a function of a in case of U = sin(27rz) + 
(l/4)sin47rz and symmetric velocity modulation w(z) as in 
Fig. [21 Other parameters are «o = 10, T = 1, and different 
5 as indicated. Current inversion at 5 = tt/2 occurs at a = 
a c = 0.05, as indicated by the vertical dashed line and a filled 
square. Inset shows U(z). 



Here we explored the possibility of rectifying their mean 
current with the help of a spatially modulated propul- 
sion speed as realizable, for example, in the experiments 



of Refs. 
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Whereas in a static potential with sym- 
metric barriers the mean current is zero, it becomes non- 
zero and can efficiently tuned with a spatially modulated 
self-propulsion and a phase shift against the potential. 
This is reminiscent to the well-known Seebeck ratchets 
[2r| of passive systems. We were able to analytically de- 
termine the mean current in the adiabatic limit of slow 
rotational diffusion for any temperature and compared 
it to numerical results. Finally we demonstrated that 
the modulated self propulsion can have a profound effect 
on the mean current in a ratchet potential. It largely 
enhances the current or can even reverse it with only a 
weak modulation. 

The rectification discussed here is universal and does 
not crucially depend on temperature. Therefore, it 
should also be found in non-Brownian active agents, such 
as the run-and-tumble particles of Ref. [23|. Our re- 
sults should be testable with the active Janus particles 
of Refs. 2(J 21 1, where the self-propulsion speed can be 
tuned by light intensity. In particular, the adiabatic limit 
of slow rotation mainly used in this article is accessible 
in the experiments. The rotational decorrelation time 
I?,? 1 due to thermal diffusion depends on the third power 
of the particle radius and can be varied in the range 
D^ 1 = 1 . . . 100s. Together with achievable velocities of 
Vo = 1/im/s and typical spatial periods of L w 30/im [2fj| . 
the adiabatic limit condition D~ Y = 100s ^> L/vq = 30s 
is easily met. This makes the system ideal for explor- 
ing the control and rectification of active particles with 
modulated self-propulsion speed. 

We thank C. Bcchingcr and V. Volpe for helpful dis- 
cussions. This work was partly supported by the re- 
search training group GRK 1558, funded by Deutsche 
Forschungsgcmeinschaft . 
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EFFECTIVE LANGEVIN DYNAMICS IN THE 3D 

CASE 

In the ovcrdampcd limit, the coordinate r — (x, y, z) 
of the particle obeys the Langevin equation [l[ 

r = -p,d z U{z)e z + v {z)p + £(t), p = 77(f) x p, (1) 

where p stands for the translational mobility, e z is the 
unit vector of the z-axis and the unit vector p points 
along the instantaneous direction of the ballistic mo- 
tion. £(f) and 77(f) represent translational and rota- 
tional noise, respectively. The later are linked to the 
absolute temperature T and the rotational diffusivity 
D r via the fluctuation-dissipation relations (£(f)£(t')) = 
2pk B TS(t - f) and (r)(t)r)(t')) = 2D r S(t - f). 

If the multiplicative noise in the second equation in 
Eqs. ([T]) is treated in the Stratonovich interpretation, the 
corresponding Smoluchowski equation for the probability 
density p(r,p,t) becomes 

% + divJ t = D r R 2 p, (2) 
at 

with the translational current Jt = —jid z U{z)pe z + 
vo(z)pp — /ifcsTVp and the rotational operator R = 
pxV p | 

It can be readily shown that 



R 2 



V p (p 2 V p - p(p • V p )) 



(3) 



Taking into account the one-dimensional nature of the 
modulation function vo{z) and of the potential U(z), and 
by using the spherical coordinates for the direction vec- 
tor p = (sin 6 cos a, sin 9 sin a, cos 9) with 9 representing 
the angle between the z-axis and the unit vector p, the 
Smoluchowsky equation Eq. ([2]) can be written in terms 
of the distribution density that only depends on one spa- 
tial coordinate p(z, 9, a, f ) 



dp 
dt 



dz 



D r d 
shifldfl 



sin 9 



dp 
39 



D r d 2 P 
sin 2 9 da 2 ' 



(4) 



with the z-component of the translational current J z = 
[—pd z U(z) + vo(z) cos9 — pkBTd z ]p(z, 9, a, t). 

The Smoluchowski Eq. ((4]) can be integrated over the 
polar angle a to yield 



dp dJ z D r d ( . dp 
dt + ~dz~ ~ ^9 89 l Sm d9 



(5) 



With the change of variables q = cos 9, the Smolu- 
chowski Eq. ([5]) is transformed to 



dp I dJz = p 9 
dt dz r dq 







(W) 

d_ 

dq 



2^dp 



D r qp + D r yjl- q 2 — sj\ - q 2 p 



(6) 



with p = p(z, q, t). 

The expression in the round brackets in the r.h.s. of 
Eq. (|6|) represents the (negative) probability current j q of 
the stochastic process q(t) 



d r 



jq = D r qp + D r \J\ - q 2 — y/1 - q 2 p 



(7) 



which corresponds to a Langevin equation in the 
Stratonovich interpretation 



q = -D r q + ^l-q 2 ri{t), 



(8) 



where 77(f) is the white Gaussian noise with the correla- 
tion function (r](t)r](t')) = 2D r S(t - f). 

The stationary distribution R(q) is found by setting 
the probability current to zero, i.e. j q = 0. This yields 



R(q) = const = -, ?G[-1,1]. 



(9) 



In order to compute the conditional average (q(t) \ qo, to), 
we use the time-dependent Fokker-Planck equation for 
the process q(t), which is obtained from the Smolu- 
chowski equation Eq. ([6]) by integrating the later over the 
coordinate z in the limits z <E (—00, 00) 



dp dj q 
dt dq 



0. 



(10) 



Thus, multiplying Eq. (|10[) with u and integrating by 
parts, we have 



d_ 

dt 



(q(t) I qo,t ) 



1 dp{q,t) j 

q — — d q 

. x dt 



jgdq, (11) 



where we have used the fact that the time-dependent 
probability current j q vanishes at q = ±1. Using Eq. ([7]) 



2 



we obtain 
d 



at 



q(t) I io,t ) 



-D r (q(t) | q , t ) 

D r J y/l-q^y/l-q^pdq 



= -2D r (q(t) | q ,t ). 
The general solution of the Eq. (|12p is given by 
(q(t)\q 0> t )=q e- 2D ^- t °\ 



(12) 



(13) 



This results can now be used to compute the station- 
ary auto-correlation function (ACF) (q(t)q(t')). It is well 
known that for a stationary Markov process, the ACF is 
given by [|| 



q(t)q(t')) = / (q(t) \ q , t')q R(q ) dq 



L e -2D r (t-t') 



(14) 



and find according to Q the stationary conditional den- 
sity p s {z\ q) and the corresponding current J(q) 



Ps(z\q) 



pk B T 



I+[L) 



\ _ e — v Lq/tik B T 



J{q) = fik B T 



I + (L)I-(L) 

— p — v a Lq/ 'fJ,k B T 



-HQ) 



where the factor 1/3 before the exponent represents the witl1 ^(z) = fo E±1 (y)dy and I(q) 



stationary variance (qfys — (1/2) /_■, <?o ^<?o = V^- 



^ E- 1 {z)I+{z)dz and 



p-V"r\ I ^Qff(z,g)-fQZg 



(16) 



FINITE-T ADIABATIC CURRENT 

In the limit of slow rotation, i.e. when q(t) changes 
slowly with time, we set q(t) = q — const in the equation 

dU eS (z,q(t)) 
z = v q(t) — hf(i), I 15 ) 
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